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Abstract 

We describe the present status of the pion distribution amplitude (DA) as it origi- 
nated from several sources: (i) a nonperturbative approach, based on QCD sum rules 
with nonlocal condensates, (ii) an 0{as) QCD analysis of the CLEO data on F^'^ ^{Q'^) 
with asymptotic and renormalon models for higher twists, and (iii) recent high-precision 
lattice QCD calculations of the second moment of the pion DA. Then we show the com- 
parison of the results for the pion electromagnetic form factor, obtained in analytic 
perturbation theory, with JLab data on Ft^{Q'^). After that we introduce the improved 
model for nonlocal condensates in QCD and show its consequences for the pion DA and 
77* TT transition form factor. In order to facilitate possible applications of BMS and 
improved "bunches" we suggest approximate analytic descriptions of their boundaries. 
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1 Generalized QCD SRs for mesonic distribution am- 
plitudes 

The pion distribution amplitude (DA) parameterizes the matrix element of the nonlocal axial 
current on the light cone pj 



{0\d{z)-f,j,E{z,0)u{0)\n{P)) 



i 



(1) 



The gauge-invariance of this DA is guarantied by the Fock-Schwinger string 



zg / A^{T)dT^ 



E{z,0) = Vexp 



, inserted in between separated quark fields. Physical meaning of this DA is quite evident: it 
is the amplitude for the transition vr(P) — *■ u{Px) + d{P{l — x)). It is convenient to represent 
the pion DA using expansion in Gegenbauer polynomials Cn'^'^{'2x — 1), which are one-loop 
eigenfunctions of ERBL kernel [2], [3] 



^.(x; ft') = v^%x) \l + Y, «2n(/i') C',i\2x - 1) 



(2) 



n>l 



where {p^^{x) = 6x{l — x) is the famous asymptotic DA. This representations means 
that all scale dependence in (y9^(a;;/i^) is transformed to the scale dependence of the set 
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{a2(/i^), a4(/i^), . . .}. We mention here that ERBL solution at the 2-loop level is also possi- 
ble with using the same representation ([2]) [U El El [7j . 

In order to construct reliable QCD SRs for the pion DA moments one needs, as has 
been shown in [HE], to take into account the nonlocality of QCD vacuum condensates. For 
an illustration of the nonlocal condensate (NLC) model we use here the minimal Gaussian 
model 

(g(0)g(^)) = (gg)eH^'l^^/« (3) 

with a single scale parameter = (A;^), which characterizes the average momentum of 
quarks in the QCD vacuum. Its value has been estimated in QCD SR approach and also on 
the lattice [ini nulla [13]: 

Aj = 0.35 - 0.55 GeV^ . (4) 

We see that A^ is of an order of the typical hadronic scale ^0.6 GeV^. 

Let us write down as an example the NLC QCD SR for the pion DA ip-jrix). To produce 
it one starts with a correlator of currents J^^ix) and Jl^.jyiO) = d{0) h'-f^ {n'V)^u{0) with 
light-like vector n, = 0, obtains SRs for the moments {x^)tj and then applies the inverse 
Mellin transform, (x^)^ =^ V^7r(a^)- As a result we obtain 

i=2V,3L,4Q 

with A = Ag/M^. The local limit A — ^ of this SR is specified by the appearance of 
5-functions concentrated at the end-points x = and x = 1, for example, (y94Q(x;A) = 
9[S{x) + 5{1 - x)]. 
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Figure 1: We show here contributions to the Eq. (EI) due to perturbative 
loop (dotted line) and due to 4Q-condensate: (p^^Q^x) — in standard QCD 
SRs, and ^p^^^ix, = 0.55 - 0.80 GeV^) — in NLC QCD SRs. 

The minimal Gaussian model ([3]) generates the contribution ip4Q{x; A) shown in Fig. [1] in 
comparison with the perturbative one for the standard (local) and the NLC types of the 
SR. We see that due to completely different behaviour of perturbative and condensate terms 
in the local QCD SR case it is difficult to obtain some kind of consistency. Alternatively, 
the NLC contribution is much more similar to the perturbative one. Just by this reason 
we have very good stability in the NLC SR case! After processing SR ([5]) for the moments 
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Figure 2: Panel (a): Moments (^^)7r with = 2, . . . , 10, obtained using the NLC 
SR ([5]), are shown by triangles with error-bars in comparison with the asymptotic 
DA moments (squares). Panel (b): The allowed values of parameters 02 and 
of the bunches (jH]), evaluated at /x^ = 1.35 GeV^ for three values of the nonlocality 
parameter = 0.4,0.5, and 0.6 GeV^. 



(C^)7r = Jq ^Trix) {2x — 1)^ dx, we restore the pion DA ip-,r{x) by demanding that it should 
reproduce these first five moments (^*)7r, z = 2 , 4 , . . . , 10 with using the minimally possible 
number of the Gegenbauer harmonics in representation ([2]). It appears that NLC SRs for 
the pion DA produce a bunch of self-consistent 2-parameter models at /ig ^ 1.35 GeV^: 



^r'^l^; /^o) = V^'"^^) [1 + «2(/i^) C!^\2x - 1) + a.ifil) Cf\2x - 1)J . (6) 

The central point corresponds to af^^ = +0.188, 04^^ = —0.130 in the case where = 0.4 
GeV^, whereas other allowed values of parameters 02 and 04 are shown in the left panel of 
Fig. [2] as the slanted rectangle [H] . By self-consistency of these solutions we understand 
that all of them produce, in accord with ([6]), the inverse moment of the pion DA 

(x-^)^""^"^ = 3.17 ±0.20. (7) 

And this range is in a good agreement with the estimation dictated by the special SR for 
this moment, which can be obtained through the basic SR by integration in x with the 




Figure 3: The inverse moment {x ^)^, obtained using the NLC SR is shown 
by the solid line (central value) with error-bars, shown as dashed lines. 
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corresponding weight x ^ (at /Xq — 1-35 GeV^): 

(x-i)^^ = 3.30 ±0.30, (8) 

see Fig. [3l It is worth to emphasize here that the moment could be determined only 

in NLC SRs because there are no end-point singularities. 

Comparing the obtained pion DA with the Chernyak&Zhitnitsky (CZ) one pTSj reveals 
that although both DAs are two-humped they are quite different — BMS DA is strongly 
end-point suppressed, as illustrated in Fig. Hl^a). To display this property more explicitly we 
show in the panel (b) of this figure the comparison of BMS and CZ contributions of different 
bins to inverse moment (a;~^),r, calculated as J^^"*""^^ (p{x)dx and normalized to 100%. 
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Figure 4: Panel (a): We compare here shapes of three DAs — BMS (solid line), 
CZ (dashed line), and the asymptotic DA (dotted line). Panel (b): Histograms 
for contributions of different bins to inverse moment (x~^),r are shown for CZ and 
BMS DAs. 



2 Analysis of CLEO data on F^^-^iQ'^) pion DA 

Many studies [HI [T71 [HI [TSj [TH [20] have been performed in the literature to determine the 
pion DA using the high-precision CLEO data [21] on the pion-photon transition form factor 
F^^*^(Q^). In particular, in [19] we have used Light-Cone Sum Rules (LCSR) [221 US] to 
the next-to-leading-order accuracy of QCD perturbation theory to examine the theoretical 
uncertainties involved in the CLEO-data analysis in order to extract more reliably the first 
two non-trivial Gegenbauer coefficients 02 and 04, which parameterize the deviation from 
the asymptotic expression (p^^. 

Why does one need to use Light-Cone SRs (LCSRs) in analyzing the experimental data 
on 7* ((5)7(5') 7r°-transition form factor? For 3> m^, -C the QCD factorization 
is valid only in the leading twist approximation and the higher twists are of importance [23] . 
The reason is quite clear: if — one needs to take into account interaction of a real photon 
at long distances of order of 0(1/ \f(f ). Then in order to account for long-distance effects in 
perturbative QCD one needs to introduce the light-cone DA of the real photon. Instead of 
doing so, Khodjamirian [22] suggested to use the LCSR approach, which effectively accounts 
for long-distances effects of the real photon using quark-hadron duality in the vector channel 
and dispersion relation in q^. Schmedding and Yakovlev realized this approach to analyze 
the CLEO data on the 7*7 it transition form factor with the next-to-leading-order (NLO) 
accuracy of perturbative QCD part of LCSR [16j . 
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We improved the NLO analysis of the CLEO data by taking into account the following 
points: (i) the NLO evolution for both v^(a;,Qgxp) (^s{Qlxp) with accurate taking into 
account heavy quark thresholds; (ii) the relation between the "nonlocality" scale and the 
twist-4 magnitude 5|^_4 ^ Ag/2 was used to re-estimate (5t„_4 = 0.19±0.02 at = 0.4 GeV^; 
(iii) the possibility to extract constraints on {x~^)t, from the CLEO data and to compare 
them with those we have from NLC QCD SRs. 

Results of our analysis [19] are displayed in Fig. [5l Solid lines in all figures enclose 
the 2(T-contours, whereas the la-contours are enclosed by dashed lines. The three slanted 
and shaded rectangles represent the constraints on (02, 04) posed by the QCD SRs for 
corresponding values of = 0.6, 0.5, 0.4 GeV^ (from left to right). All values are evaluated 
at /i^ = 2.4 GeV^ after the NLO evolution. 




(1') «i 



Figure 5: Three 2a- and la-contours (solid and dashed lines, correspondingly) of the 
admissible regions following from the analysis of the CLEO data for different values 
of 6^: (a) - for A^ = 0.4 GeV^ and 6^ = (0.19 ± 0.02) GeV^; (b) - for A^ = 0.5 GeV^ 
and 52 = (0.235 ±0.025) GeV^; (c) - for A^ = 0.6 GeV^ and 5^ = (0.29 ±0.03) GeV^ 

We see that the CLEO data definitely prefer the value of the QCD nonlocality parameter 
Ag = 0.4 GeV^. We also see in Fig. [5]^c) (and this conclusion was confirmed even with 20% 
uncertainty in twist-4 magnitude, see also Fig. [6]^a)) that CZ DA (■) is excluded at least 
at 4a-level, whereas the asymptotic DA (♦) — at 3a-level. In the same time our DA (X) 
and most of the bunch (the slanted green-shaded rectangle around the symbol H) are inside 
la-domain. Instanton-based Bochum (ik) and Dubna (A) models are near 3a-boundary and 
only the Krakow model pi], denoted in Fig. [6]^a) by symbol 4-, is close to 2a-boundary. 




Figure 6: (a): The results of the CLEO data analysis for the pion DA parameters 
((x-^)^"P/3 - 1, evaluated at /i^ ^ 1 GeVl (b): LCSR predictions for Q^F^*^^^{Q^) 
for the CZ DA (upper dashed line), BMS- "bunch" (shaded strip), and the asymptotic 
DA (lower dashed line) in comparison with the CELLO (diamonds, [25]) and the 
CLEO (triangles, [21]) experimental data, evaluated with the twist-4 parameter value 
52 = 0.19 GeV^ [IS] and at fily = 5.76 GeVl 
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In Fig. M^a) we demonstrate the la-, 2a- and 3cr-contours (solid, dotted and dashed 
contours around the best-fit point (+)), which have been obtained for values of the twist-4 
scale parameter 5tw-4 ~ [0.15 — 0.23] GeV^. As one sees from the blue dashed line within the 
hatched band, corresponding in this figure to the mean value of {x~^)^/3 — 1 and its error 
bars, the nonlocal QCD sum-rules result with its error bars appears to be in good agreement 
with the CLEO-constraints on {x~^)'^^ at the Icr-level, Moreover, the estimate {x^^)^ is 
close to (x^^)™/3 — 1 = 0.24 ± 0.16, obtained in the data analysis of the electromagnetic 
pion form factor within the framework of a different LCSR method in [26l |27]. These three 
independent estimates are in good agreement to each other, giving firm support that the 
CLEO data processing, on one hand, and the theoretical calculations, on the other, are 
mutually consistent. 

Another possibility, suggested in j28], to obtain constraints on the pion DA in the LCSR 
analysis of the CLEO data - to use for the twist-4 contribution renormalon-based model, 
relating it then to parameters 02 and 04 of the pion DA. Using this method we obtain [29] 
the renormalon-based constraints for the parameters a2 and 04, shown Fig. [8] in the form of 
Icr-ellipse (dashed contour). 



3 Dijet E791 data, pion form factor and CEBAF data 



Our findings are further confirmed by the E791 data [30] on diffractive dijet 7r-|-A-production. 
This is illustrated in Fig. El^a). The main conclusion here is that all considered pion DAs 
are consistent with the data, with tiny preference to the BMS DA. Indeed, following the 
convolution procedure of [31] we found [19] the following values of for three types of pion 
DAs: 12.56 (asymptotic DA), 10.96 (BMS bunch), and 14.15 (CZ DA). 
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Figure 7: (a): Comparison with the E791 data on diffractive dijet production of the 
BMS "bunch" (shaded strip), the asymptotic DA (solid line), and the CZ (dashed 
line) model, using the convolution approach of [31] • (b): The scaled pion form factor 
calculated with the BMS bunch (shaded strip) and asymptotic DA (dashed lines) 
including nonperturbative uncertainties from NLC QCD SRs [2] and renormalization 
scheme and scale ambiguities at the 0(a^)-level [32]. The experimental data are taken 
from [33] (diamonds) and [31], [35] (triangles). 

It is worth to mention here also the results of our analysis of the pion electromagnetic 
form factor using NLC dictated pion DA and Analytic Perturbative QCD [32]. These results 
are in excellent agreement with CEBAF data on pion form factor, as illustrated in the Fig. 
[7](b), where the green strip includes both the NLC QCD SRs uncertainties, generated by our 
bunch of the allowed pion DA, and by the scale-setting ambiguities at the NLO level. 
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From the phenomenological point of view, the most interesting resuh here is that the 
BMS pion DA [Tl] (out of a "bunch" of similar doubly-peaked endpoint-suppressed pion 
DAs) yields to predictions for the electromagnetic form factor very close to those obtained 
with the asymptotic pion DA. Conversely, we see that a small deviation of the prediction for 
the pion form factor from that obtained with the asymptotic pion DA does not necessarily 
imply that the underlying pion DA has to be close to the asymptotic profile. Much more 
important is the behavior of the pion DA in the endpoint region x — > , 1. 

4 New lattice data and pion DA 

Rather recently new high-precision lattice measurements of the pion DA second moment 
= Jq (2x — dx appeared [Ml ETj. Both groups extracted from their respec- 

tive simulations, values of 02 at the Schmedding-Yakovlev scale /Xgy around 0.24, but with 
different error bars. 

It is remarkable that these lattice results are in striking agreement with the estimates of 
02 both from NLC QCD SRs fT4] and also from the CLEO-data analyses — based on LCSR — 
[inilin], as illustrated in Fig. [Ht^a), where the lattice results of [37j are shown in the form of a 
vertical strip, containing the central value with associated errors. Noteworthily, the value of 
02 of the displayed lattice measurements (middle hne of the strip) is very close to the CLEO 
best fit in (+). 

5 Improved model for NLCs and pion DA 

The quark-gluon-antiquark condensates are usually parameterized in the following form 

{Q{0h,,{-9^Ay))q{x)) = {y^Xy- g^^{yx))Mi{x^,y'^,{y - xf) 

+ {Vt^Vu - 9^uy^)M2{x^, l/^ {y - x)^) , 
{q{0h5lf,{-9A^{y))q{x)) = iSf^uy^Msix^, y^, {y - x)'^) , 

with (Ai,2,3 = ^(-i,2,|)) 

00 

M.{x',y',z') = A, jjjdad(3d^ U{a,(3,^)e^'^^"^^y"^^'")'\ (9) 


The minimal model of nonlocal QCD vacuum suggests the following Ansatze 

/, (a, /5, 7) = 5 (a - A) 5 (/3 - A) 5 (7 - A) (10) 

with A = |Ag and faces problems with QCD equations of motion and gauge invariance of 
2-point correlator of vector currents. In order to fulfil QCD equations of motion exactly and 
minimize non-transversity of — correlator we suggest [38j the improved model of QCD 
vacuum with 

(«, /?, 7) = (1 + X,d, + Y,dy + YA) S {a - xA) 6 {f3 - yA) 6 (7 - zA) , (11) 

where 2; = y, A = and 

Xi = +0.082; X2 = -1.298; X3 = +1.775; a; = 0.788; (12a) 
Fi = -2.243; 1^2 = -0.239; Kj = -3.166; ?/ = 0.212. (12b) 
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Then the NLC sum rules [5] produce [38] a ^^bunch" of 2-parameter pion DA models ([6]) at 
/i^ = 1.35 GeV^. The coordinates of the central point 4- are 02 = 0.268 and 04 = —0.186. 
These values correspond to (^x~^)^'^^'^^ = 3.24 ± 0.20, which is in agreement with the result 
of an independent sum rule, viz., = 3.40 ± 0.34. In order to facilitate possible 

applications of our "bunches" we suggest the following approximate analytic descriptions of 
their boundaries — we specify upper and lower curves, 04 (02), and left and right vertical 



lines, fixed by their abscissa values, Og' : 

a4^^'+(a2) = 0.177 -2.41 02 + 7.77 a^- 14.0 a^; (13a) 

a^^^'~ia2) = -0.027 - 0.246 02- 4.12 0^ + 8.86 a^; (13b) 

BMS, L niQ/i BMS, R noKi /to \ 

02 = 0.134; ttg = 0.251, (13c) 

for the BMS bunch, whereas for the improved bunch 

a^P'+(a2) = 0.352 -4.25 02 + 14.0 a^- 18.9 a^; (14a) 

a^P'-(a2) = 0.588 - 6.89 a2 + 20.7 a^- 23.6 a^; (14b) 

a^P'L = 0.177; a^^'^ = 0.346. (14c) 



Allowed values of both bunch parameters 02 and 04 after NLO-evolution to /x^ = 5.76 GeV^ 
are shown in Fig. [8] in a form of shaded slanted rectangles around central points H and 
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Figure 8: The results of the CLEO data analysis for the pion DA parameters 02 and 
04, evaluated at /igy = 5.76 GeV^. The lattice results of [37] are shown for comparison 
as shaded area, whereas the renormalon-based la-ellipse of [29] is displayed by the 
green dashed line. The shaded strip in panel (a) shows the corresponding results 
for the BMS- "bunch", whereas in panel (b) — the new bunch, corresponding to the 
improved Gaussian model of QCD vacuum. 

We emphasize here that BMS model [14], shown in Fig. [8](a) by symbol J<, is inside 
the allowed region dictated by the improved QCD vacuum model. This means that all the 
characteristic features of the BMS bunch are also valid for the improved bunch: NLC-dictated 
models are end-point suppressed, although are double-humped. 

We see in Fig. Mjo) that the improved bunch [38] is even in a better agreement with the 
recent lattice results [37], shown in the form of a vertical strip, containing the central value 
with associated errors. Remarkably, the value of 02 of the displayed lattice measurements 
(middle line of the strip) is very close to the central point of the "improved bunch" ("♦"), 
whereas the whole bunch, dictated by the improved NLC QCD SRs [38], is inside the stripn 
and completely inside the standard CLEO la-ellipse. 

^The same is valid for the BMS "bunch" as well, as can be seen in Fig. [5Ja). 
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6 Conclusions 

Let me conclude with the following observations: 

• NLC QCD SR method for the pion DA gives us the admissible bunches of DAs for 
each value of Xq. 

• NLO LCSR method produces new constraints on the pion DA parameters (02 and a^) 
in conjunction with the CLEO data. 

• Comparing results of the NLC SRs with new CLEO constraints allows to fix the value 
of QCD vacuum nonlocality: = 0.4 GeV^. 

• This bunch of pion DAs agrees well E791 data on diffractive dijet vr + A-production, 
with JLab F(pi) data on the pion form factor, and with recent lattice data. 

• Taking into account QCD equations of motion for NLCs and transversity of vacuum 
polarization allows us to put the pion DA bunch just inside Icr-ellipse of CLEO-data 
constraints. 
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